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Abstract 

We develop a systematic method for computing a renormalized light-front field theory 
Hamiltonian that can lead to bound states that rapidly converge in an expansion in free- 
particle Fock-space sectors. To accomplish this without dropping any Fock states from the 
theory, and to regulate the Hamiltonian, we suppress the matrix elements of the Hamiltonian 
between free-particle Fock-space states that differ in free mass by more than a cutoff. The 
cutoff violates a number of physical principles of the theory, and thus the Hamiltonian is 
not just the canonical Hamiltonian with masses and couplings redefined by renormalization. 
Instead, the Hamiltonian must be allowed to contain all operators that are consistent with 
the unviolated physical principles of the theory. We show that if we require the Hamiltonian 
to produce cutoff-independent physical quantities and we require it to respect the unviolated 
physical principles of the theory, then its matrix elements are uniquely determined in terms of 
the fundamental parameters of the theory. This method is designed to be applied to QCD, but, 
for simplicity, we illustrate our method by computing and analyzing second- and third-order 
matrix elements of the Hamiltonian in massless <fi 3 theory in six dimensions. 
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1 Introduction 



It may be possible for hadron states to rapidly converge in an expansion in free-particle Fock-space sectors 
(free sectors) in Hamiltonian light-front quantum chromodynamics (HLFQCD). This will happen if the 
Hamiltonian satisfies three conditions. First, the diagonal matrix elements of the Hamiltonian in the basis 
of free-particle Fock-space states (free states) must be dominated by the free part of the Hamiltonian. 
Second, the off-diagonal matrix elements of the Hamiltonian in this basis must quickly decrease as the 
difference of the free masses of the states increases. If the Hamiltonian satisfies these first two conditions, 
then each of its eigenstates will be dominated by free-state components with free masses that are close to 
the mass of the eigenstate. The third condition on the Hamiltonian is that the free mass of a free state 
must quickly increase as the number of particles in the state increases. If the Hamiltonian satisfies all 
three conditions^], then the number of particles in a free-state component that dominates an eigenstate 
will be limited from above by the constraint that the free mass of the free-state component must be close 
to the mass of the eigenstate. This means that the eigenstate will rapidly converge in an expansion in 
free sectors^. 

We argue that if one uses Hamiltonian light-front field theory (HLFFT) rather than equal-time field 
theory, then it is possible to calculate a Hamiltonian that satisfies these conditions. If the field theory 
is to be complete, then the Hamiltonian must satisfy these conditions without a truncation of the Fock 
space or the removal of particle-number-changing interactions. Based on the early work of Dyson [1] 
and Wilson [2] , and the more recent work of Wegner [3] and Glazek and Wilson [4] , there have been a 
substantial number of efforts to derive such Hamiltonians perturbatively [5-12]. However, a calculation 
of the QCD Hamiltonian beyond second order in perturbation theory requires one to take the scale 
dependence of the coupling into account, and none of these earlier methods demonstrates how to do this. 
In this paper we develop an alternative method for calculating a renormalized Hamiltonian that satisfies 
the conditions discussed above. We do not truncate the Fock space or remove particle-number-changing 
interactions. We show how the scale dependence of the coupling is taken into account in general and we 
present explicit examples that demonstrate how this is done in practice. 

If we are to be able to derive a Hamiltonian that satisfies the conditions discussed above, then we must 
work in an approach in which it is possible for the vacuum to be dominated by few-body free sectors. 
This does not seem possible in equal-time field theory unless the volume of space is severely limited. In 
HLFFT, however, each particle in the vacuum has a longitudinal momentum of exactly zero[]. This means 
that we can force the vacuum to be empty by requiring every particle to have a positive longitudinal 
momentum. The effects of the excluded particles must be replaced by interactions in the Hamiltonian, 
and a calculation of these interactions may have to be nonperturbative. We view the calculation of these 
interactions as an extension of our method and plan to include them in the future. Until we do consider 

1 Exactly how quickly the Hamiltonian's off-diagonal matrix elements must decrease and the mass of a free state must 
increase is not known; so we assume that the rates that we are able to achieve are sufficient. This must be verified by 
diagonalizing the Hamiltonian. 

2 The coefficients of the expansion for highly excited eigenstates may grow for a number of free sectors and then peak 
before diminishing and becoming rapidly convergent. 

3 This is because there are no negative longitudinal momenta and momentum conservation requires the three- momenta 
of the constituents of the vacuum to sum to zero. 
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these interactions, we are avoiding the vacuum problem rather than solving it. 

In light-front field theory, all the dynamics of the Hamiltonian are in the invariant-mass operator and 
their matrix elements are trivially related. For convenience, we work with the invariant-mass operator 
directly and refer to it as the Hamiltonian. To see how we can derive a Hamiltonian that satisfies our 
conditions, we divide it into a free part and an interacting part: 

M 2 (A) = M) + Mj(A). (1) 

The interacting part of the Hamiltonian is a function of the regulator, A. The free states are eigenstates 
of the free Hamiltonian: 

M)\F)=M%\F). (2) 

The diagonal matrix elements of the Hamiltonian are given by 

(F\M 2 (A)\F) = (F\M}\F) + {F\Mj(A)\F) 

= M 2 F (F\F) + {F\Mj(A)\F). (3) 

We assume that the Hamiltonian can be computed perturbatively, in which case the free part of the 
Hamiltonian will dominate the diagonal matrix elements. This fulfills the first condition on the Hamilto- 
nian. 

We regulate the Hamiltonian with a cutoff A such that 

-(Ml, -M?) 2 

{F\M 2 (A)\I) ~e ^ , (4) 

and then as \Mp — M 2 \ grows, the off-diagonal matrix elements of the Hamiltonian quickly diminish, 
thus fulfilling the second condition on the Hamiltonian. We show in Appendix A that in HLFQCD it is 
reasonable to expect that Mj, grows as least quadratically with the number of particles in \F), perhaps 
much faster. This fulfills the third and final condition on the Hamiltonian. 

The cutoff violates a number of physical principles of the theory, such as Lorentz invariance and 
gauge invariance. This implies that the Hamiltonian cannot be just the canonical Hamiltonian with 
masses and couplings redefined by rcnormalization. Instead, the Hamiltonian must be allowed to contain 
all operators that are consistent with the unviolated physical principles of the theory. The key point of 
this paper is that if we require the Hamiltonian to produce cutoff-independent physical quantities and 
we require it to respect the unviolated physical principles of the theory, then its matrix elements are 
uniquely determined in terms of the fundamental parameters of the theory. Although the focus of this 
paper is the computation of the Hamiltonian, any observable can be calculated in our approach. 

Our main assumption is that in an asymptotically free theory the Hamiltonian can be determined 
perturbatively, provided the cutoff is chosen so that the couplings are sufficiently small. Although we com- 
pute the Hamiltonian perturbatively, it can be used to obtain nonperturbative quantities. For example, 
bound-state masses can be computed by diagonalizing the Hamiltonian. The disadvantage of pertur- 
bative renormalization is that nonperturbative physical quantities will be somewhat cutoff-dependent. 
The strength of this cutoff dependence has to be checked after nonperturbative physical quantities are 
computed and is not considered here. 
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In general, field theories have an infinite number of degrees of freedom. However, since our Hamiltonian 
will cause hadron states to rapidly converge in an expansion in free sectors, approximate computations 
of physical quantities will require only finite-body matrix elements of the Hamiltonian. In addition, 
since we assume that we can renormalize the Hamiltonian perturbatively, we do not implement particle 
creation and annihilation nonperturbatively in the renormalization process. This allows us to work with 
a finite number of degrees of freedom when calculating the required finite-body matrix elements of the 
Hamiltonian. 

We are going to apply the method to QCD in future work; however, the complexities of QCD are 
largely irrelevant to its development. For this reason, we choose to illustrate the method with a much 
simpler theory, massless (f> 3 theory in six dimensions. This theory has some similarities to QCD that make 
it a suitable testing ground. In particular, its diagrams have a similar structure to QCD diagrams and 
it is asymptotically free. Nonperturbatively the theory is unstable; however, we only use it to illustrate 
how to compute the Hamiltonian perturbatively, and so the instability is irrelevant. 

In Section 2 we specify how the Hamiltonian is regulated and restrict it to produce cutoff-independent 
physical quantities. In Section 3 we restrict the Hamiltonian to respect the physical principles of the 
theory that are not violated by the cutoff. In Section 4 we show how the restrictions that we place 
on the Hamiltonian allow us to uniquely determine its matrix elements in terms of the fundamental 
parameters of the theory. Section 5 contains example calculations of a few second-order matrix elements 
of the Hamiltonian. Section 6 contains a perturbative example of how the cut-off Hamiltonian fulfills 
the restriction that it must lead to cutoff-independent physical quantities. In Section 7 we calculate the 
matrix element of the Hamiltonian for (j> — > <jxf> through third orderQ and demonstrate how the coupling 
runs in our approach. Finally, in Section 8 we provide a summary and a short discussion of the extension 
of our method to QCD. 

2 Regulation and Renormalization 
2.1 The Cutoff 

We work in the basis of free states. We write the Hamiltonian as the sum of the canonical free part and 
an interaction (see Appendix B for conventions and definitions): 

M\h) = M 2 f + Mj(A). (5) 

At this point the interaction is undefined. To force the off-diagonal matrix elements of the Hamiltonian 
to decrease quickly for increasingly different free masses, we regulate the Hamiltonian by suppressing its 
matrix elements between states that differ in free mass by more than a cutoff: 

(F\M 2 (A)\I) = {F\M)\I) + {F\M 2 {k)\I) 

= M 2 F {F\I)+e-— {F\V{K)\I), (6) 

4 We neglect contributions to this matrix element that depend on the cutoff only through the Gaussian regulating factor. 
See Sections 4 and 7 for details. 
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where \F) and \I) are eigenstates of the free Hamiltonian with eigenvalues Mp and Mf, and A FJ is the 
difference of these eigenvalues: 

A FI = M 2 - M 2 . (7) 

V(A) is the interaction with the Gaussian factor removed, and we refer to it as the "reduced interaction." 
To determine the Hamiltonian, we must determine the reduced interaction. 

Assuming (F\ V(A) \I) does not grow exponentially as A FI gets large, the exponential in Eq. (||) 
suppresses each off-diagonal matrix element of the Hamiltonian for which A FI is large compared to A . 
This regulates the Hamiltonian and forces it to weakly couple free states with significantly different free 
masses. 



2.2 The Restriction to Produce Cutoff-Independent Physical Quantities 



Unfortunately, any regulator in HLFFT breaks Lorentz invariance, and in gauge theories it also breaks 
gauge invariance^]. This means that in HLFFT renormalization cannot be performed simply by redefining 
the canonical couplings and masses. Instead, the Hamiltonian must be allowed to contain all operators 
that are consistent with the unviolated physical principles of the theory. In addition, since there is 
no locality in the longitudinal direction in HLFFTQ, these operators can contain arbitrary functions of 
longitudinal momenta. 

To uniquely determine the operators that the Hamiltonian can contain, as well as the coefficients 
of these operators, we place a number of restrictions on the Hamiltonian. The first restriction is that 
the Hamiltonian has to produce cutoff-independent physical quantities. We impose this restriction by 
requiring the Hamiltonian to satisfy 



M 2 {A) = U(A, cA) M 2 (cA) U f (A, cA), 



(8) 



where c > 1 and is otherwise arbitrary, and U is a unitary transformation that changes the Hamiltonian's 
cutoff. Note that we are considering M 2 to be a function of its argument; i.e. A4 2 (cA) has the same 
functional dependence on cA that M 2 (A) has on A. To see that Eq. (^) implies that the Hamiltonian will 
produce cutoff-independent physical quantities, note that the Hamiltonian on the right-hand side (RHS) 
of the equation can be replaced by iterating the equation: 

M 2 (A) = U(A,cA)M 2 (cA)U f (A,cA) 

= U{A, cA) U(cA, c 2 A)j M 2 (c 2 A) [[/ f (cA, c 2 A) £7* (A, cA) 

= ~U(A, cA) U(cA, c 2 A) U(c 2 A, c 3 A)j M 2 (c 3 A) \u\c 2 A, c 3 A) C/+(cA, c 2 A) C/t(A, cA) 

(9) 

5 Our regulator breaks these symmetries because the mass of a free state is neither gauge-invariant nor rotationally 
invariant (except for transverse rotations). 

6 That there is no longitudinal locality in HLFFT is evident from the fact that the longitudinal momentum of a free 
particle appears in the denominator of its dispersion relation, (p? + m 2 )/p+. 
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Thus, since c > 1, Eq. (g) implies that M 2 (A) is unitarily equivalent to liniA^oo _M 2 (A). This means 
A4 2 (A) will produce cutoff- independent physical quantities^. 
For simplicity, we define 

A' = cA, (10) 

and then Eq. (||) takes the form 

M 2 (A) = U(A,A') M 2 (A') f/t(A, A'). (11) 

To calculate the Hamiltonian perturbatively, we need to be able to implement this equation pertur- 
batively. In our approach, the cutoff that regulates a Hamiltonian also serves to define the scale of the 
coupling that appears in the Hamiltonian; so the coupling in A4 2 (A) is defined at the scale A. To use Eq. 
( pi] ) perturbatively, we have to perturbatively relate the coupling at the scale A to the coupling at the 
scale A'. This perturbative relationship is well-defined only if A' is not very large compared to A [13]. To 
fulfill this requirement and Eq. ( |l0| ) (recall c > 1), we choose A' to satisfy 

A < A' < 2A. (12) 

Now Eq. (|ll]), which restricts the Hamiltonian to produce cutoff- independent physical quantities, can be 
imposed perturbatively. 

2.3 The Unitary Transformation 



To make Eq. (|llj) a complete statement, we have to define the unitary transformation. The transforma- 
tion is designed to alter the cutoff implemented in Eq. (|^) , and is a simplified version of a transformation 
introduced by Wegner [3], modified for implementation with the invariant-mass operator. The transfor- 
mation is uniquely defined by a linear first-order differential equation: 

dU(A,A') 



d(A-*) 



T(A)U(A,A% (13) 



with one boundary condition: 



f/(A,A) = l. (14) 

It is shown in Appendix C that U(A, A') is unitary as long as T(A) is anti-Hermitian. We define 

T(A)=[M^ 2 (A)], (15) 

which is anti-Hermitian. The transformation introduced by Wegner uses 

T(A)= [Ml{A) 1 M 2 {A)} 1 (16) 

where 7Wq(A) is the part of M 2 (A) that is diagonal in the free basis. (It includes the diagonal parts of all 
interactions.) Our transformation regulates changes in eigenvalues of -M- 2 , and Wegner 's transformation 
regulates changes in eigenvalues of M^A). 

7 If the theory is not asymptotically free, then this argument cannot be used in perturbation theory. 



G 



To solve for A4 2 (A) perturbatively, we need to turn Eq. (|ll| ) into a perturbative restriction on the 
reduced interaction, V(A). We begin by taking the derivative of Eq. (pi]): 

= [[M},M 2 (A)],M 2 (A)], (17) 

where Eq. ([l3]) and its adjoint and Eq. Jl5| ) are usedQ. Since the free part of the Hamiltonian commutes 
with itself and is independent of the cutoff, 

?0W = [[M},M 2 I (A)] ) M 2 (A)]. (18) 

We assume that M 2 (A) can be written as a convergent (or at least asymptotic) expansion in powers 
of Mf(A'): 

OO 

M 2 (A) = Y. via) > ( 19 ) 

a=0 

where O ([A^j(A')] a ), i/ ' — M 2 , and we must determine the higher-order i/^'s. Substituting 

Eq. ( [Tg| ) into Eq. ([l8]), matching powers of AA 2 {A') on both sides, and doing some algebra yields 

^ = i[[M> {b) }^ a - b ^ (20) 

where the sum is zero if a = 0. This equation tells us how the O ( [.Mj(A')] °) contribution to M. 2 (A) 
depends on the cutoff in terms of lower-order contributions, and can be used to iteratively calculate 
A"1 2 (A) in powers of M 2 (A'). 

Taking a matrix element of Eq. (^p| ) between free states \i) and \j) for the case a = 1, we obtain 

4 )A ?r (21) 



d(A- 4 ) 13 13 
Solving this equation with the boundary condition 

M 2 {A)\ =M\A') (22) 



leads to 



(i) 

v>j = exp 



I A— A' 
1 1 

A 71 ~ A 4 



Mj(A%. (23) 



Mj(A') is regulated such that M. 2 (A')ij is proportional to exp(— Af^-A' -4 ), which means that Eq. ( f23|) 
shows that vQ' is proportional to exp(— A^- A~ 4 ). 

The higher-order ?/ a )'s can be found the same way that was found, by taking matrix elements 
of Eq. ( |20| ) and using the lower-order v^'s and the boundary condition, Eq. (p2|). When products of 
Mj(A') are encountered, complete sets of free states have to be inserted between adjacent factors. The 
«( a )'s that result from this procedure are all proportional to exp(— A^A~ 4 ), which shows that U{A, A') 
does indeed take a Hamiltonian with a cutoff A' and produce one with a cutoff A. 

If one computes v {2) and v (3 \ then Eq. @ gives M 2 {A) in terms of Mj(A') to O ([M 2 (A')]' 
Using this and 

8 The restrictions that we place on the Hamiltonian in the next section prohibit A4 2 (A') from depending on A. 
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(F\Mj(A)\I)=e-^ (F\V(A)\I), (24) 

which follows from Eq. (^), one can show that the perturbative version of Eq. ([TIJ) in terms of the 
reduced interaction is 

V(A)-V(A') = SV, (25) 
where SV is the change in the reduced interaction and is a function of both A and A': 
{F\SV\I) = ^WVMWWVMWT^&K' 1 ) 

K 

+ \Y, (F\ V(A') \K) (K\ V(A') \L) (L\ V(A') \I) Tt' A '\F, K, L, I) 

K.L 

+ o([V(A')f). (26) 
In this equation, the sums are over complete sets of free states and the cutoff functions are defined by 

T^' A '\F,K,I) = f-L- - -M />A'- 4 A PX A Ki _ e 2A^A™A K; j (2?) 
\ A FK A KlJ ^ ' 



and 



T^ AA) (F,K,L,I) 



A FK A Kh) \ A LI A FL 



(e 2A '~ 

2A'~ 4 A FJf A KL ^ e 2A- 4 A Fi A LJ _ ^2A'~ 4 A FL A LJ ^ 

(e 2A '~ 



A KL A LI J \ A KI A FK 

A FK + A IK ( , e 2A'- 4 (A FK A if/ +A KI/ A LI ) _ g 2A- ; , .A , h A KI +A KI A, 

*KL ^LI / A KL A LI + A FK A KI 



A FL + A IL f e 2A'-*(A FK A KL +A FL A LI ) _ g 2A 1 i A , K A + A , , A. ,, 

A FK A KLJ A FK A KL + A FL A LI 



(28) 



The above definitions for the cutoff functions assume that none of the A's that appear in the denominators 
is zero. In the event one of them is zero, the appropriate cutoff function is defined by 

I\ (A ' A,) (A = 0) = lim T; (A,A,) (A). (29) 

To summarize, Eq. ( |25| ) is the perturbative version of the restriction that J\A 2 (A) has to produce 
cutoff- independent physical quantities, expressed in terms of the reduced interaction and the change in 
the reduced interaction. 



3 Restrictions on the Hamiltonian from Physical Principles 

Eq. ( p5| ) is the first restriction on the Hamiltonian. To uniquely determine the Hamiltonian, we need to 
place additional restrictions on it, and we do this using the physical principles of the theory that are not 
violated by the cutoff. 
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3.1 Symmetry Properties 

Any LFFT should exhibit momentum conservation, boost invariance, and transverse rotational invariance. 
Our cutoff does not violate any of these principles; so we restrict the Hamiltonian to conserve momentum 
and to be invariant under boosts and transverse rotations. 



3.2 Cluster Decomposition 



Cluster decomposition is a physical principle of LFFT that is partially violated by our cutoff. However, 
we can identify the source of the violation and still use the principle to restrict the Hamiltonian. 

To see how to do this, note that momentum conservation implies that any matrix element of the 
reduced interaction can be written as a sum of terms, with each term containing a unique product of 
momentum-conserving delta functions. For example, (0304| V(A) \4>i4>2) can be written in the form 

647r 5 (^ 



(0 3 4 | V(A)\<t>lfo) = 



+ 
+ 



Pi)S (5) {pi +P2 -P3 -Pi)\ F^( Pl , P2 , P3 , P4 ,A) 

6Att 5 p+S^{ Pi -p 3 )] [6Att 5 p+5^{ P2 -p 4 )] F& (p u p 2 , p 3 , p±, A) 
6A7r 5 p+6 (5 \ Pl - Pi )} \m^p+8^\ P2 -p 3 )l F^( Pl ,p 2 ,p 3 , P4 ,A). 



(30) 



We have included a factor of 647r 5 and a longitudinal momentum factor with each delta function because 
our normalization of states produces these factors naturally. Cluster decomposition implies that the fW's 
in Eq. ( |30| ) cannot contain delta functions of momenta [15]. It also implies that each term in Eq. ( |30| ) 
has exactly one delta function associated with the conservation of momenta of interacting particles, with 
any additional delta functions associated with the conservation of momenta of spectators. This means 
that the first term on the RHS of Eq. ([30j) is the contribution from four interacting particles, and the 
second term on the RHS can have two contributions: one for which 0i and 03 are spectators, and one 
for which <p 2 and 04 are spectators. Similarly, the third term on the RHS can have two contributions: 
one for which 0i and 04 are spectators, and one for which 02 and 03 are spectators. 

Normally, contributions to matrix elements depend on spectators only through momentum-conserving 
delta functions and the corresponding factors of longitudinal momenta, in which case we can write 

P2-P3~Pd F (1) (p 1 ,p 2 ,p 3 ,p4,A) 

+ [647r 5 p+S^( Pl - P3 )\ [6iTT 5 p+S^( P2 - P4 )\ [F^ 1 \p 2 ,p i ,K) + F^ 2 \ Pl ,p 3 ,K) 

+ [e^p+S^ipt-pi)] [u^ P +5^{p 2 -p 3 )\ [f ( - 3 < 1 \ P2 , P3 ,A) + F^ 2 H Pi , P4 ,A)~\, (31) 

where we explicitly show the contributions from different sets of spectators. However, we use a cutoff 
on differences of free masses of states, and the free mass of a state does not separate into independent 
contributions from each particle in the state. A particle's contribution to the free mass of a state is 
f 2 /x, where x is the fraction of the state's longitudinal momentum carried by the particle and rj_ 
is the transverse momentum of the particle in the state's center-of-mass frame. When the longitudinal 
momentum of any other particle in the state changes, x changes because the total longitudinal momentum 
of the state changes. Thus our cutoff partially violates the cluster decomposition principle, such that the 
can depend on P + , the total longitudinal momentum of the states: 



<&04| V(A)|<M2> = Mir 5 ( P j 



pt)S {5 Hpi 
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(030 4 | V(A) |0! 2 ) = 

64 7 r 5 p+'5 (5) bi -Ps) 
647r 5 p+<5 (5) (pi -p 4 ) 



64^ 5 (p+ +p+)^ 5 )( Pl +p 2 - pa -p 4 )] F^( Pl ,p 2 ,p 3 , Pi , A) 
64^+^ 5 )(p 2 - p 4 )] [f( 2 ' 1 )(p2, P4 , P + , A) + F^fpi.ps, P + , A) 
64^ 5 p+^ 5 ) (pa - p 3 )] [F^ (pa , p 3 , P+ , A) + F< 3 ' 2 > (pi , p 4j P+ , A) 



■(32) 



The result of this analysis can be generalized to formulate a restriction on the Hamiltonian: when a 
matrix element of the reduced interaction is written as an expansion in the possible combinations of 
momentum-conserving delta functions, the coefficient of any term in the expansion depends only on the 
momenta of the particles that are interacting in the term, the total longitudinal momentum of the states, 
and the cutoff. 

It is worth mentioning that if one uses a cutoff on free masses of states rather than differences of free 
masses, then the F^'s in Eq. © can be functions of all the momenta in the matrix element, which 
is a signal of a much more severe violation of cluster decomposition. If one uses a cutoff on free energy 
differences, then cluster decomposition is not violated at all, but longitudinal boost invariance is lost. 



3.3 Transverse Locality 

Ideally, a LFFT Hamiltonian should be local in the transverse directions, and thus each of its matrix ele- 
ments should be expressible as a finite series of powers of transverse momenta with expansion coefficients 
that are functions of longitudinal moment aQ. In our case, the cutoff suppresses interactions that have 
large transverse momentum transfers and replaces them with interactions that have smaller transverse 
momentum transfers. This is equivalent to suppressing interactions that occur over small transverse 
separations and replacing them with interactions that occur over larger transverse separations; so we do 
not expect our interactions to be perfectly transverse-local. Nonetheless, we expect that interactions in 
.M 2 (A) should appear local relative to transverse separations larger than A -1 or, equivalently, to trans- 
verse momenta less than A. This means that for transverse momenta less than A we should be able to 
approximate each matrix element of .M 2 (A) as a finite power series in p±/A. We enforce this by assuming 
that transverse locality is violated in the weakest manner possible, i.e. that any matrix element of the 
Hamiltonian can be expressed as an infinite series of powers of transverse momenta with an infinite radius 
of convergence. 

3.4 Representation of the Theory of Interest 

In the remainder of this section, we place a number of restrictions on the Hamiltonian that limit it to 
describing the particular LFFT of interest, massless </> 3 theory. We work in six dimensions so that the 
theory is asymptotically free. 

We assume that we can compute the Hamiltonian perturbatively, which means that we can expand 
V"(A) in powers of the coupling at the scale A. Our cutoff has no effect in the noninteracting limit; so 

9 This series is actually multiplied by a product of momentum-conserving delta functions which has no impact on the 
locality of the Hamiltonian. 
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our Hamiltonian must reproduce free massless scalar field theory in this limit. According to Eq. (|6j), this 
means that V(A) vanishes in the noninteracting limit. 

In massless (ft 3 theory, the only fundamental parameter is the coupling; so we require the Hamiltonian 
to depend only on it and the scale. In this case, the expansion of V(A) takes the form 



f(a) = £>;fc-)(a), 



(33) 



where g A is the coupling at the scale A. We refer to V^(A) as the 0{g r K ) reduced interaction, although 
for convenience the coupling is factored out. 

g A is the correct fundamental parameter for massless (ft 3 theory if and only if its definition is consistent 
with the canonical definition of the coupling. The canonical definition is 

-l 



64tt 5 p+5 (5) (pi -P2 



P3) 



(34) 



We choose 



64Tr 5 P fd (5 \ Pl -ft -ps)l {fc(ft 3 \M 2 (A) |^) 



64ir 5 p+S (5 \ Pl - pa - p 3 j\ (4> 2 <ft a \V(A) |<M 



I P2± =P3±=0 ; pf =pj = 3P+ 



(35) 



Our definition of the coupling is consistent with the canonical definition because the conditions on the 
matrix elements in Eq. (|3^) have no effect on the matrix element in Eq. (0), which is momentum- 
independent. According to Eq. (|33"|), the Hamiltonian is coupling coherent [14] because the couplings 
of its noncanonical operators are functions only of the fundamental parameters of the theory and they 
vanish in the noninteracting limit. 

If the Hamiltonian is to produce the correct second-order scattering amplitudes for massless (ft 3 theory, 
then the first-order reduced interaction must be the canonical interaction: 

= (2n) 5 P+ J D 1 D 2 D 3 [ala 1 a 2 5^\p 3 ~p 1 -p 2 ) + alala 1 S^(p 2 +p 3 -pi)]- (36) 

A proof of this statement exists, but the proof is tedious, and so we neglect to present it. 

A perturbative scattering amplitude in (ft 3 theory depends only on odd powers of the coupling if the 
number of particles changes by an odd number and only on even powers of the coupling if the number 
of particles changes by an even number. We require our Hamiltonian to produce perturbative scattering 
amplitudes with this feature. 

In the remaining sections, we show how the restrictions that we have placed on the Hamiltonian 
uniquely determine it in terms of the fundamental parameters of the theory and allow it to produce 
correct physical quantities. As a check on our procedure, we can verify that the Hamiltonian that 
results from the procedure satisfies the restrictions that we have placed on it. The restrictions that 
we use are based on a subset of the physical principles of the theory. The remaining principles, such as 
Lorentz invariance and gauge invariance (in a gauge theory), must be automatically respected by physical 
quantities derived from our Hamiltonian, at least perturbatively. If they are not, then they contradict 
the principles we use and no consistent theory can be built upon the complete set of principles. 
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4 The Method for Computing Matrix Elements of the Hamil- 
tonian 



In this section we develop a general formalism for using the restrictions that we have placed on the 
Hamiltonian to calculate its matrix elements. To begin, we consider the restriction that forces the 
Hamiltonian to produce cutoff-independent physical quantities: 

V(A)-V(A') = SV. (37) 

This restriction is in terms of the reduced interaction and the change in the reduced interaction. 

SV is defined in Eq. (^6|), which makes it clear that since V(A') can be expanded in powers of g , , so 
can SV: 

oo 

SV = J29l Sv(t) - (38) 

t—2 

We refer to SV^ as the 0(g t r ) change in the reduced interaction, although for convenience the coupling 
is factored out. Note that SV^ is a function of A and A'. 
Now Eq. ( |37|) can be expanded in powers of g A and g A , : 

oo oo oo 

E jy® ( a) - e 1 y {t) ( a') = E & 6V(t) ■ (39) 

t=l t=l t=2 

This equation is a bit tricky to use because it involves the coupling at two different scales. To see how 
they are related, consider the matrix element of Eq. ( |37| ) for ifa — > </>2</ > 3 : 

(020 3 | V(A) |0!> - (020 3 | V(A') |0!> = {<h(fo\6V \<h) • (40) 

According to the definition of the coupling, this equation implies 

9 A -9 A , = [64ttV<5 (5) ( Pl - P2 - P3 )\ ^ (fofo \8V\4n) I , n + + t + - (41) 

Since V^ 1 ) changes particle number by 1, inspection of Eq. ( p6|) reveals that (0203 1 ^ |</>i) is 0(g^,); so 

g A =g Al +0(g 3 A ,). (42) 

This implies 

oo 

■9a =.? v +E^^( A ' A '). ( 43 ) 

s=3 

where the C s 's are functions of A and A'. In Section 7 we calculate C3 explicitly. For an integer t > 1, 
Eq. (|43|) implies 

00 

■9l-.i+E^ + ^( A ' A '), ( 44 ) 

where the -B tjS 's are also functions of A and A', and can be calculated in terms of the C s 's by raising Eq. 
(ph to the t th power. 
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We substitute Eq. (Q into Eq. ( p9|) and demand that it hold order-by-order in g A ,. At 0{g r l ) 
(r > 1), this implies 

r-1 

F W (A) - ^ (r) (A') = SV {r) - J2 B r - s , s V {r - s) (A), (45) 

where SV^ = 0, and we define any sum to be zero if its upper limit is less than its lower limit. 

We have restricted the Hamiltonian to respect approximate transverse locality, which means that 
its matrix elements can be expanded in powers of transverse momenta. This means that the matrix 
elements of (A) can also be expanded in powers of transverse momenta. Each term in a transverse- 
momentum expansion of a matrix element {F\ (A) \I) is either cutoff-dependent or cutoff-independent. 
We define Vqq(A) to be the cutoff-dependent part of V^(A), i.e. the part that produces the cutoff- 
dependent terms in transverse- momentum expansions of matrix elements of V^(A). We define to 
be the cutoff- independent part of V^ r )(A), i.e. the part that produces the cutoff- independent terms in 
transverse- momentum expansions of matrix elements of V^(A). Then 

7W(A) = ^(A)+7(?. (46) 

When we substitute this equation into Eq. (|45|), the cutoff-independent parts of the terms on the left-hand 
side (LHS) cancel, leaving 

r-1 

V&>(A) - V^(A') = 5VV - £ B r ^M r - s \A). (47) 

s=2 

This equation can be used to calculate the O(sl) reduced interaction in terms of the lower-order reduced 
interactions. Note that for r = 1, the RHS of the equation is zero, implying that the cutoff-dependent 
part of is zero. Inspection of Eq. ( |36| ) shows this to be the case. 

In the remainder of this section, we summarize the results of Appendix D, which contains more details 
and rigor than are necessary here. 

Recall that momentum conservation implies that any matrix element V(A) \I) can be written as 
an expansion in unique products of momentum-conserving delta functions. This means that an arbitrary 
matrix element of Eq. (|47]) can be expanded in products of delta functions and thus is equivalent to a set 
of equations, one for each possible product of delta functions. Given approximate transverse locality, each 
of the resulting equations can be expanded in powers of transverse momenta. Matching the coefficients of 
powers of transverse momenta on either side of these equations allows us to rigorously derive the following 
results (see Appendix D for details). 

First, the cutoff-dependent part of the O(g^) reduced interaction is given in terms of lower-order 
reduced interactions by 

r-1 

(F\5V^\I)-Y,Br- s , s (F\V( r - s HA)\I) , (48) 

s= 2 J A terms 



(F\V^(A)\I) = 



where "A terms" means that the RHS is to be expanded in powers of transverse momenta and only the 
terms in the expansion that depend on A are to be kept. Recall that SV is defined in Eq. (|26|). 

Second, the cutoff-independent part of V^IA) is the part with three interacting particles and no 
transverse momentum dependence. If there is no such part, then V^ r )(A) is completely determined by 
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Eq. ©. Third, V^(A) can have a cutoff-independent part only if r is odd. Fourth, the coupling runs 



at odd orders; i.e. C s is zero if s is even [see Eq. (43)]. Fifth, there is no wave function renormalization 
at any order in perturbation theory in our approach because this would violate the restrictions that we 
have placed on the Hamiltonian. 

The sixth and final result from Appendix D is that the cutoff- independent part of the O(g^) re- 
duced interaction for odd r > 3 is given in terms of the cutoff-dependent part and lower-order reduced 
interactions by the integral equation^] 

1 



(F\V£\I) = 



-Br, 2 



{F\5V {r+2) \I) 



Ext. k 



r+1 

E 

s=3 



(F\V& +2 -^\I) 



(49) 



- 3 int. part. 



"Ext. k± — > 0" means the limit in which the transverse momenta in the external states are taken to zero. 
"3 int. part." means the products of momentum-conserving delta functions that appear on the RHS are 
to be examined and only the contributions from three interacting particles are to be kept. 

Eq. (|4^) is an integral equation because Vqj is nested with lower-order reduced interactions in 
integrals in 8V^ r+2 \ A cursory examination of the definition of 5I/( r+2 ) may lead one to believe that 
Eq. ( p9[ ) is useless because it requires one to know V^ r+l \A'). However, since r is odd, V^ r+1 ^ (A') = 
V^ r D +1) (A'), and the dependence of 5V^ r+ ^ on V^ r D +1) (A') can be replaced with further dependence on 



Vq{ and Vq'^(A') (and lower-order reduced interactions) using Eq. ( |4S| ) with r — * 

In the remaining sections, we apply the results of this section to calculate a number of example 
matrix elements of the Hamiltonian, demonstrating that they are indeed uniquely determined in terms 
of the fundamental parameters of the theory by the restrictions that the Hamiltonian must produce 
cutoff-independent physical quantities and be consistent with the unviolated physical principles of the 
theory. 



r{r) 



1. 



5 Second- Order Matrix Elements 

In this section we illustrate our method by calculating some second-order matrix elements of the Hamilto- 
nian. The matrix elements of the Hamiltonian are defined in terms of the matrix elements of the reduced 
interaction by Eq. (Jfy. To calculate matrix elements of the Hamiltonian to second-order, we need to cal- 
culate the corresponding matrix elements of V ^ (A) , the second-order reduced interaction. As mentioned 
in the previous section, for V"' r ^(A) to have a cutoff-independent part, r must be odd. According to Eq. 
(Es), this means that V^(A) is given in terms of the second-order change to the reduced interaction by 



(F\V {2 \A)\I) = (F\SV^\I) 



(50) 

A terms 



where the second-order change to the reduced interaction is defined in terms of the first-order reduced 
interaction and the cutoff function: 

{F\SV^\I} = -^(F\V (1) \K)(K\V {1) \I)T { 2 ^ A '\F,K,I). (51) 

K 



10 It is very difficult to prove that integral equations of this type have a unique solution; so we simply assume it is true in 
this case. 
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Figure 1: A diagrammatic representation of (02 1 SV^ |0i}, a matrix element of the second-order change 
to the reduced interaction. The numbers label the particles. 

5.1 Example: ((p 2 \M 2 (A) |0i) 

As a first example, we compute the matrix element (02 1 A^ 2 (A) \<fii) to O(g^). To compute this matrix 
element, we use \I) — |0i) and \F) — 1 2 ) ■ Then the matrix element of the second-order change to the 
reduced interaction is 

(0 2 |^ (2) |<^i) = ^(0 2 |FW|ir)(^|F( 1 )|0 1 )Ti A ^(F,ir,J). (52) 

1 K 

Since changes particle number by 1, \K) has to be a two-particle state and Eq. ( f32| ) becomes 

(021^10!) = \ [d 3 D 4 (0 2 |F« 10304} (030 4 | V^\^)T^ h '\F,KJl (53) 



where the completeness relation in Eq. (177) is used (see Appendix B for conventions and definitions) 
and \K) = |0304). Here (0 2 |<5V( 2 ) |0i) is represented diagrammatically in Fig. 1. Using the expression 
for in Eq. @, 

(0 2 |<ST/ (2) |0i} = - A J ' D 3 D 4 S 2M S ltU T^' A '\F,K,I). (54) 

The integral over p 3 can be done with one of the delta functions. At this point it is convenient to change 
variables from p^ and p 4 ± to x and r± : 

Pi = (pt ,P4_l) = {xpf,xp 1± +f±), (55) 

where x is the fraction of the total longitudinal momentum carried by particle 4 and r± is the transverse 
momentum of particle 4 in the center-of-mass frame. Since p^ is constrained, we do not display it in the 
list of components of p 4 . The momentum-conserving delta functions imply 

P3 = ([1 -x]Pi, [1 -x]pi± -fx) ■ (56) 

Note that all longitudinal momenta are positive, although we do not explicitly show these limits in 
integrals. This implies, for example, that < x < 1. With the change of variables, Eq. ( |54]) becomes 

(fa\SVW\<h) = \piS^(p 2 - Pl )Jd i r A _dx-^— ) T^ A '\F,K,I). (57) 

We are using massless particles; so the free masses of the states \I) = |0i) and \F) = |02) are zero: 

M] = M 2 F = 0. (58) 

The particles in the intermediate state \K) — 10304} are also massless, but the state still has a free mass 
from the particles' relative motion: 
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Then from the definition of the cutoff function T% A ' A \ 



Tt' A '\F,K,I) = - 2<1 ~ X) (e- 2 ^^^ _ e - 2A - 4 ^fe^j . (60) 

After the integral is done, Eq. (|57| ) becomes 

WVW\^) = Ki^LssJl (A 2 -A' 2 ). (61) 
Applying Eq. ( |50| ) yields the matrix element of the second-order reduced interaction: 

(0 2 |^C 2 )(A)|0i) = A^A 2 ^!-^. (62) 

Using the definition of the Hamiltonian in terms of the reduced interaction, we find that to second-order, 
the matrix element of the Hamiltonian is given by 

{h \M 2 (A)\cb 1 )^S 1 , 2 A" 9 ^ YE L- 3] ^. (63) 

This matrix element of the Hamiltonian has been completely determined to second-order, in terms 
of the fundamental parameters of the theory, by the restrictions that the Hamiltonian has to produce 
cutoff-indcpcndcnt physical quantities and has to respect the unviolated physical principles of the theory. 
If we wanted to compute physical quantities, this is one of the matrix elements of the Hamiltonian that 
we might need. We would have to choose A and determine g A by fitting data since the theory contains 
one adjustable parameter. Note that Eq. ( ]63| ) is consistent with all the restrictions that we have placed 
on the Hamiltonian. 

5.2 Example: (0 3 4 | M 2 (A) \(f>i<fr 2 ) 

Next we calculate the matrix element {0 3 c^>4 1 A4 2 (A) |0i02) to O(g^). We use \I) = |0i02) and |F) = 
10304). Then the matrix clement of the second-order change to the reduced interaction is 

(fafaldvW |0i0 2 ) = i^(0304|U«|A)(A-|T/W|^ 1 ^ 2 )T 2 (A ' A ' ) (F,^,/). (64) 

Z K 

Since changes particle number by 1, \K) has to be either a one- or three-particle state. When 
V' 1 ' is substituted into the RHS of this equation and all the creation and annihilation operators from 
the possible intermediate states and the interactions are contracted, we see that {fo^SVM |0i0 2 ) has 
contributions from a number of different processes. We define (0304 1 SV^ 10102)^ to be the contribution 
to (0 3 4 | 5V^ 10102) from the process shown in Fig. 2-i, where 1 < i < 9. Then 

9 

(030 4 | SV& |0i0 2 ) = (^304| SV (2) |0i0 2 ), . (65) 
i=l 

Eq. (^) suggests that we define a contribution to the matrix element of the second-order reduced 
interaction for each contribution to (0304 1 SV^ |0i02): 
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(2-1) 








(2-8) 



(2-9) 



Figure 2: A diagrammatic representation of the contributions to (03^4 1 SV^ \4>i<j>2), a matrix element of 
the second-order change to the reduced interaction. The numbers label the particles. 
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(0304| K (2) (A) 1^102), = (foMSV^^ifo), 



A terms 



and then 



/>3^4| K (2) (A) |01<^2> = ^3^1 ^ (2) (A) 10102), 



(66) 



(67) 



For each contribution to (0304 1 5V^ 2 ' |0i02), we use the change of variables 

Pi = {xP + ,xP±_ +fj_), 

P2 = ([1 - x]P+, [1 - x]P± - rl), 

P3 = (yP + ,yP ± + q±), 

p A = ([l-y]P+,[l-y}P ± -q ± ), 



(68) 



where P = p± + pi . Then the initial and final free masses are 



M 2 F 



x(l — x) 



(69) 



5.2.1 The Annihilation Contribution 

The contribution to (0304 1 SV^ |0i02) from the annihilation process, shown in Fig. 2-1, is 

(030 4 |<5^ 2 )|0 1 2 ) 1 = 5/ ^5 (0304| W 1 )|0 5 )(0 5 |^ 1 )|0 1 2 )T 2 (A ' A ' ) (F, if,/) 



1/1 



2.34 



2 V Mj, Mf 



-2A'~ 4 MfMf 



-2A" 4 MpM, 2 



(70) 



where the intermediate state is |if) = |05). According to Eq. (|66|), the corresponding part of the matrix 
element of the second-order reduced interaction is 



(0 3 4 |^ 2 )(A)|0 1 2 ) 1 



Si 



1/1 



2,34 



1 



2 VMf, Mf 



-2A' ~ 4 Mj,M] 



-2A~ 4 Mj.Aff 



(71) 



A terms 



If we expand everything multiplying the delta function in powers of transverse momenta, the lowest 
order terms from the two exponentials cancel and leave two types of terms: those depending on A and 
those depending on A'. We can remove the A'-dependent terms without altering the cancellation of the 
lowest term by replacing the first exponential with a 1: 



(0 3 4 |^ (2) (A)|0 1 2 ) 1 



1 / 1 



J 12,34 



2 V Mf ' Mj 



1 - e 



-2A~ 4 Mj,Mf 



(72) 
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5.2.2 The Two-Particle Self-Energy Contributions 



The rest of the contributions to (0304| SV^ |0i02) have a three-body intermediate state, \K) — |050607). 
The contribution from the two-particle sclf-cncrgy, shown in Fig. 2-2, is given by 



(0304|<5V (2) |0i0 2 ) 2 - - A — I D 5 D 6 D 7 T 2 (A ' A ' } (F, K, I) 8 3 ,m S 7A 8i fiB 5 7 , 2 - 

Axy J 

We use the change of variables 

p 5 = (zpf,zpi± + k±), 
and do the integrals over p 6 and p 7 using the delta functions. Then we find that 

Pe = ([l-z]pt,[l-z]pix-k±), 
and the free mass of the intermediate state is 



kl 



x(l — x) z(l — z)x ' 



(73) 



(74) 



(75) 



(76) 



and 



6 3 4 | 5V^ |0102) 2 = n B1 _R ^1,3 #2,4 



2 64?r 5 

<5l,3 ^2,4 



d k\dz- 



-4 A 2 

ik _ e 



-2A""A 



1536tt 



(77) 



According to Eq. (|66J), the corresponding part of the matrix element of the second-order reduced inter- 
action is 

1 HF 
1536tt 3 V 2' 



(0304 | ^ (2) (A)|0 1 



2)2 



$1,3 $2,4 A 2 



(78) 



Note that Eqs. (62) and ( f78[ ) indicate that to 0(g 2 ) in this theory the only effect a spectator has on the 
self-energy is to produce an extra delta function. This is what one would expect if cluster decomposition 
were maintained. However, there is no guarantee that this will hold to all orders, and it is known that 
in gauge theories (0304 1 V^ 2 ^(A) |0i02) 2 depends on / P + . This is allowed because our cutoff partially 
violates cluster decomposition. 



According to Fig. 2, 

<030 4 |V^(A) 10x02)3 

(0 3 4 |^ 2 )(A)|0 1 2 ) 4 
(0 3 4 |^ 2 )(A)|0 1 2 ) 5 



(030 4 |F (2) (A)|0 1 2 ) 2 

(030 4 |^ (2) (A)|0102) 2 
(030 4 |^ (2) (A)|0102) 4 



3^4 
l<->2 

3^4 
3^4 



Sl,3 S 2A A 2 



1 



1536tt 3 V 2 : 
1 



1536tt 3 V 2 : 
1 



1536tt 3 



(79) 
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5.2.3 The Exchange Contributions 



The exchange contribution to the matrix element of the second-order change in the reduced interaction, 
shown in Fig. 2-6, is 



(0 3 4 |^ 2 > \<f> 1( f, 2 ) 6 = 1 — ! _ J D 5 D 6 D 7 T^ A '\F,K,I)S it67 5 7 , 2 §i,m h,z 

1 1 



2,34 



2 x - y 

where the free mass of the intermediate state is 



Ml 



^KF ^KI 



-2A' 4 A KF A KI _ e -2A i A KF A KI 



(80) 



(81) 



1 — x y x — y 

Eq. ( |66| ) implies that the corresponding part of the matrix element of the second-order reduced interaction 
is given by 



(0304|^ (2) (A)|010 2 ) 6 = -<5l2,34 



1 1 



1 



1 



^KF ^KI 



1 - e 



-2A" 4 A,^ c ,A l 



(82) 



2x-y 

According to Fig. 2, the remaining parts of the matrix element of the second-order reduced interaction 
are given by 

= (0304| V (2) (A)|0 1 2 ) £ 



<0 3 4 |^ (2) (A) |010 2 ) 7 

(030 4 | V (2 >(A) |0r 2 ) g 

(0 3 4 |^ (2) (A)|0 1 2 ) 9 



(030 4 | y (2) (A) |0i0a>« 

(0304|^ (2) (A)|010 2 ) 6 

(0 3 4 |^ (2) (A)|0i0 2 ) 8 



3—4 
1—2 

3—4 
3—4 



(030 4 |V^(A)|0 1 2 ) £ 
(030 4 | y (2) (A)|0!0 2 ) £ 



v— >l— y , q±^>-q± 

x — >1 — x . rj_ — > — r_L 

y— l— a , q±^-q± 



(83) 



£ — >l — x , f± — > — r_L 



5.2.4 The Complete Matrix Element and Transverse Locality 



Using Eqs. (|6j) and (33), we can write the full second-order matrix element of the Hamiltonian: 



6 3 4 | M 2 (A) |010 2 



[Si. 3 S 2A + S 1A «5 2 ,3)Ml + e- A_4A " 5 2 ^(0304|V r ^(A)|0 1 2 ) i . (84) 



i=l 



By inspection, it is clear that this result respects all the restrictions we have placed on the Hamiltonian, 
except perhaps the requirement that it can be represented as a power series in transverse momenta with 
an infinite radius of convergence. To see how to check this, consider the first term in the sum on the 
RHS: 



-A- 4 Ai. r „2 



'gl (030 4 | W 2 )(A)|0 1 2 ) 1 = W^- A "( M -- M? ) 2 ( 1 ^ 2 - + 1 ^ 



Mp Mj 



x 1 - e 



-2A~ 4 mIm? 



(85) 



The presence of the mass denominators might seem likely to cause a problem with the convergence of a 
power series expansion. Everything multiplying the delta function can be expressed as a power series in 
transverse momenta if it can be expressed as a power series in M F and Mf, since it is a function only of 
these and since Mp and Mf are themselves power series in transverse momenta. We can write 
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(3-1) (3-2) (3-3) 

Figure 3: A diagrammatic representation of the contributions to (020304| SV^ |0i), a matrix element of 
the second-order change to the reduced interaction. The numbers label the particles. 



1 , 1 



1 

A 1 " 1 3 



{MpMf + M F M?) - (MpMf + MpMj 



(86) 



If we approximate the RHS of this equation by keeping only N terms, the difference between this finite 
sum and the RHS of Eq. (p3) can be made arbitrarily small, for any values of Mp and Mf, by making 
N large enough. This means that e~ A Apl 9 2 (0304| V^(A) 0i02)i can indeed be represented as a 
power series in transverse momenta with an infinite radius of convergence. The other contributions to 
(0304 A / 1 2 (A) |0i02) can be analyzed similarly, and the conclusion is the same for each. 

5.3 Example: (0 2 030 4 |^ 2 (A) |0i) 

The final second-order matrix element of the Hamiltonian that we calculate is (020304 1 A4 2 (A) |0i). 
To compute this matrix element, we use \I) = |0i) and \F) — |0 2 0304). The matrix element of the 
second-order change to the reduced interaction is 

(0 2 030 4 | ^10!) = i^(0 2 3 4 |FW| J ft:)(ir|FW|0 1 )T 2 (A ' A ' ) (F,ir,/). (87) 

K 

Since V^ 1 ' changes particle number by 1, \K) has to be a two-particle state. 

We define (0 2 0304| SV^ to be the contribution to (020304| 5V^ |0i) from the process shown 
in Fig. 3-i, where 1 < i < 3. Then 

3 

(0203041 SV^ |0i) = J2 (0203041 ^ (2) \<h)i ■ 



Eq. (50) suggests that we define a contribution to the matrix element of the second-order reduced 
interaction for each contribution to (020304| SV^ |0i): 



(0203041 (A) |0i) 



(0203041 8V^ |0lV 



(89) 



A terms 



and then 



(020304| V {2) (A) |0i) = J2 (020304| V {2) (A) |0i) 4 



(90) 
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For each contribution to (020304| |0i), we use \K) — |0506) and the change of variables 

P2 = {xpt,xpi±+f±), 

P3 = {ypt,ypi± + Q±), 

Pa = {zpt,zPi± + w±), 

where 

x + y + z = 1 , 
f±+q± + w± = 0. 

Then the relevant differences of free masses are 

A, 



(91) 



(92) 



*FK — 



X FI 







WJ_ 


X 


y 


x + y 






+ .) 


\X + y 


rj_ 






X 


y 


z 



(93) 



\A pK ^KlJ ^ ' 



and the contribution to the matrix element of the second-order change to the reduced interaction shown 
in Fig. 3-1 is 

z x + y \^ FK l->. ki / 

and the corresponding contribution to the matrix element of the second-order reduced interaction is 

(020304| W 2) (A) 101)! = S h234 l— ' ' 1 N 

V K^FK ^KI J 

(95) 



X (e 2A '~ 4A FK&Ki _ e 2A- 4 A FK A JCJ ^j 



2 x + y \ A FK A KI/ 
i \ 

A terms 

As before, if we expand everything multiplying the delta function in powers of transverse momenta, 
the lowest order terms from the two exponentials cancel and leave two types of terms: those depending on 
A and those depending on A'. We can remove the A'-dependent terms without altering the cancellation 
of the lowest term by replacing the first exponential with a 1: 



(0203041^)10^ = ^34^(^-^(1 



„2A- 4 A FK -A C 



')• 



2x+y\A FK A KI JV ~ I' (96) 

According to Fig. 3, the other parts of the matrix element of the second-order reduced interaction 
are given by 



(0203041 F (2) (A)|0l} 5 



/» 2 030 4 |^ (2) (A)|0 1 ) 1 



(0203041 V {2) (A) |0!)! 



3^4 



, (97) 

y«2 , q±-<->w± 



and 



(020304| V"( 2 )(A) |0i)g - (0203041 (A) |0!)! 

= (0203041 (A) |0i)! 



2^3 , 3^4 , 4^2 
x^y , y^z , z^x 

r±— >q± , 3j_— >«?_l , wj_— >rj_ 



(98) 
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Using the definition of the Hamiltonian in terms of the reduced interaction, we can write the full second- 
order matrix element of the Hamiltonian: 

3 

(0203041 M 2 (A) |</>i) = e - A ~ 4 <«^<0 2 030 4 |V (2) (A)|0i> i . (99) 

i—l 

This result is consistent with the restrictions we have placed on the Hamiltonian. 

6 The Removal of Cutoff Dependence from Physical Quantities 

One of the restrictions on the Hamiltonian is that it has to produce cutoff-independent physical quan- 
tities. To see how this requirement is fulfilled, consider as an example the scattering cross section for 
0102 - ► 0304- At second order, the T matrix has contributions from the s, t, and u channels. The dif- 
ferent channels are linearly independent functions of the external momenta; so each contribution should 
individually have all the properties required of the full T matrix. For simplicity, we limit ourselves to 
consideration of the s (annihilation) channel. 

In our approach, the scattering matrix S is defined by 

(F\ S | J) = (F\I) - 2^5(P F - Pf) (F\ T{P T ) | J) , (100) 

where the T matrix is defined by [16] 

T(p~) = Hj(A) + Hj(A) _ 1 T(p-), (101) 
p — h + it 

and -ffj(A) is the interacting part of the Hamiltonian cut off by A: 

*m = (io2) 

To second order, 

(0 3 4 |T(pr+P 2 ~) 10102) = (0 3 04|g/(A) + g J (A) — 1 ff / (A)|0 1 2 ). (103) 

p l +p 2 - h + te 

If we neglect all noncanonical operators, the first term does not contribute and the second term can be 
calculated by inserting a complete set of free states. This leads to 

„2 

5iP F - Pf) (0304^(^7) |0 1( /. 2 ) = ^e- 2A ~ 4s2 64^ 6 )( Pl +p 2 -pa -p 4 ), (104) 
where \I) = |0i0 2 ), \F) — |0304), and s is the invariant mass-squared of |0i02). The cross section is 



cutoff-independent only if the T matrix is cutoff- independent. Thus Eq. (104) shows that the cutoff 
appears in physical quantities if noncanonical operators are neglected. If we were to use an infinite 
cutoff, then the 0(<7 2 ) contribution to the cross section would be cutoff- independent but higher-order 
contributions would be infinite. 

To test whether our procedure for computing tV( 2 (A) leads to cutoff-independent physical quantities, 
we now include noncanonical operators. This means that we have to include the first term in Eq. ( |103| ): 



6 3 4 | F/(A) |0i0 2 ) = I+^+e~ A ~ 4A "3 A 2 (030 4 | V^\K) |0 1 2 ) 1 

1 )MnH^\p 1 +p 2 -p 3 -p i ), (105) 



Pl +P2 

9 1{\ < ^ - 
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where the result for the annihilation contribution to the second-order reduced interaction is used, and 
we use the fact that we only need the on-energy-shell part. This contribution shifts the T matrix so that 
the total T matrix is 

5{P F - Pf) (0304| T(Pf) |0i0 2 ) = ^-6Att 5 S^( Pi +P2-P3-P4)- (106) 

Since g 2 is cutoff-independent at second-order [dg^/dA = O(gf)], this is the correct cutoff-independent 
result to 0(gl). 

Our restrictions on the Hamiltonian have led to its matrix elements being such that they compensate 
for the presence of the cutoff in physical quantities. No proof exists, but we expect that our procedure 
leads to a renormalized Hamiltonian that produces exactly correct perturbative scattering amplitudes 
order-by-order. When this Hamiltonian is used nonperturbatively, however, there should be some cutoff 
dependence, as we have mentioned in previous sections. 



hH M 2 {A) \(f)i) to Third Order 



We wish to calculate (0203 M 2 (A) 0i) to third order to demonstrate how a higher-order calculation 
proceeds and to derive the running coupling in our approach. In Appendix D, we deduce that matrix 
elements of even-order reduced interactions cannot have contributions from three interacting particles. 
This means that (0203 Ad 2 (A) |0i) has only odd-order contributions. The first-order contribution is 
defined in terms of the first-order reduced interaction, which we have defined in Eq. (36). The third- 
order contribution is defined in terms of the third-order reduced interaction. In this section, we compute 
the third-order contribution, neglecting the cutoff-independent part of the third-order reduced interaction. 
Eq. ( ff9|) shows that a calculation of the cutoff-independent part would be a fifth-order calculation, which 
we prefer to avoid at this time. 

According to Eq. (Eg) , the cutoff-dependent part of the third-order reduced interaction is given by 



(0 2 3 |^(A)|0 1 ) 
i 2 03| SV^ |0i) is defined by 



(0 2 3 | 8V&> |0i) - Bi, 2 (0 2 3 | |0r) 



A terms 



(107) 



<020 3 | ^ (3) |0l) 



^{^\y {1) \K) (i,|^ 2 )(A')|0i)T 2 (A ' A,) 

K 

i£^3l^ (2) (A')|A0 (K\vV\<h)T^ 

K 

\ £ (020 3 | W \ K ) (K\ VW \L) (L\ V« |0t) 7f' A '\ 

K.L 



(108) 



where we are suppressing the dependence of the cutoff functions on the states. Note that B\ 2 = C3 and 
is given by 



Bi, 



C 3 



(0 2 03 1 (H^ 3 ) |0O 



(02031 |0!)1 (020 3 |^ 3 )|0i) 



(109) 



=p 3 ±=Q wi=vi=hpi 
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Eq. (107) then becomes 



(•MalV^AJI&H 



(^3|^ (3) |^i)-(^3|^ (3) |^i) 



(110) 



A tern 



This indicates that (0203| ^cd(A) |0i) can be computed solely in terms of (0203 1 <5T^^ 3 ^ |0i), which we 
now calculate. 

All the matrix elements of the second-order reduced interaction that can appear on the RHS of Eq. 



(108) were calculated in the previous section. When (0203| SV^ |0i) is expanded in terms of the possible 
intermediate states and possible contributions to matrix elements of and V"( 2 )(A'), we see that it has 
contributions from a number of different processes. We define (0203| SV^ \(f>i) i to be the contribution 
to (0 2 3 | SV^ \4>i) from the process shown in Fig. 4-i, where 1 < i < 15. In Fig. 4, the plain vertices 
represent matrix elements of and the boxes represent contributions to matrix elements of V^(A') 
corresponding to the processes depicted in the boxes. For example, the box in Fig. 4-12 represents 
(02 ^3 1 (A') 04 ) 6 [see Fig. 2-6 and Eq. (§^)]. Note that not every contribution to every matrix 



element of V^ 2 )(A') that appears in Eq. (108) appears in Fig 4. For example, (02031 V^(A') |040s) 9 [see 
Fig. 2-9 and Eq. (|83|)] does not appear in Fig. 4, because when it is paired with (040s| T^ 1 ) |0i) and 
we integrate over and ps, it gives an identical contribution to (020 3 | V^(A') |040s} 6 . For this reason, 
(020 3 | SV^ |0i) 12 , shown in Fig. 4-12, includes contributions from both. Then (0203| SV^ |0i) is 



(0 2 3 | |0 X ) = ]T <020 3 | SV^ 100, . 

1=1 

For each contribution to (0203| SV^ 3 ' |0i), we use the change of variables 



(111) 



P2 = (ypt,ypi± + q±), 

P3 = ([i - y]pf, [1 - y]pi± - 



(112) 



which leads to the initial and final free masses 



M? 



0, 



Mp = 



y(i -y) 

7.1 The One-Particle Self-Energy Contributions 



(113) 



The contributions to the matrix element of the third-order change to the reduced interaction from one- 
particle self-energies, shown in Figs. 4-1, 4-2, and 4-3, can be combined because they have similar 
structures. For each of them we use 

\L) = |040 5 ), 
1*0 = I'M, 

(114) 



and the change of variables 
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(4-1) 



(4-2) 



(4-3) 




(4-13) (4-14) (4-15) 



Figure 4: A diagrammatic representation of the contributions to (02 03 SV^ \<fii). The numbers label 
the particles, the bare vertices represent matrix elements of and the boxes represent contributions 
to matrix elements of V^(A') corresponding to the processes depicted in the boxes. 



2G 



P5 = ([l-xlpl^l-xlfix-fx), 
which leads to the intermediate free masses 



(115) 



Mr = 0, 
M 2 L 



(116) 



x(l — x) ' 

and pe = p\. The contribution from Fig. 4-1 is 

k^MSV^Wlx = \ / <0 2 3 |^ (1) l^e) (0e|^ (1) |040s) (0405|^(^ |0x) T 3 (A ' A,) (,F, /) 



where 



^J^ dx ^T^\F,K,L,I), 



N=^pt6 {5) {pi-P2-P?,). 



(117) 



(118) 



The contribution from Fig. 4-2 is 

(0 2 3 |^ (3) |0i) 2 - \ jD 6 ( ( j> 2< p 3 \V^\ ( l> 6 )(4> e \V^\A')\ ( l> 1 }Tt' A ' ) (F,K,I) 

and the contribution from Fig. 4-3 is 

(fofo\6V^ =1 f Di D 5 <0 2 3 |W 2 >(A') 1040s), (0 4 5 | \^)Ti A ' K '\F,L,I) 



(119) 



N 
2^2 



d r±dx- 



• 4 
1 



x(l - x) \M[ M- 



l-e- 2A '- iM Z M2 A4 A ' A '\F,L,I). 



(120) 



It is easiest to combine these three contributions before doing the integral in Eq. (117). This leads to 



3 N f°° f 

]T<0 2 3 |<^ (3) |<M 4 = -J d{Ml)\e 



1 M| 2 
Mj + M| + Mj 



+ e 2K'- A M 2 L M% I 1 



Ml 



A FL Ml M 2 F A FL 



^)]-[A'-A]}. (121) 



This one-dimensional integral cannot be done analytically for arbitrary Ml; so we consider it to be a 
function of Ml that must be computed numerically. 



7.2 The Two-Particle Self-Energy Contributions 

The contributions to the matrix element of the third-order change to the reduced interaction from two- 
particle self-energies, shown in Figs. 4-4, 4-5, and 4-6, can be combined. For each of these we use 
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\L) = |04<A 5 ), 
\K) = 1060708), 



and the change of variables 

p 6 = {xp^,xp 2 ±+r±), 

P7 = ([1 - x]p£, [1 - x]p 2 ± - fj_), 

which leads to the intermediate free masses 

r2 _ 9i 



M£ = 



M 2 q 



2/(1 -y)' 

■2 



K a;(l-a;)j/' 
and 

= P2, 

Pb = Ps = Pa- 
Then these contributions are given by 
(020 3 | ^ (3) |0i) 4 = ^ J d*r ± dx x{i * x)y2 Tt' A '\F,K,L,I), 

(0203|^)„ 1>5 = ^y^^^^.^)^.^^)^^^ 

(020 3 | ^ (3) |0i) 6 ^ ^y|A' 2 Tf' A ')(F,L,I). 
Writing their sum as an integral leads to 

6 jy roc 



-2A' _4 A 2 



A, 



- e- 2A A - w - ( - ~ 7rv - , A/w 

3 





2 




" Mf 




-[A' 


M|)_ 



^(020 3 |^ 3 )|0 1 ) 4 



i=l 



2^3 



From Fig. 4, 

9 6 

^ (0 2 3 | SV& = ^(0203 1 ^ (3) |0l)i 

i=7 i=4 

but since X)i=4 (0203| SV^ |0i) i is a function only of Mj. and iV, which are invariant under 2^3 

9 6 
]T(020 3 |^ (3) |0l), - ^ (020 3 | ^ (3) |0l)i- 



j=7 i=4 
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7.3 The Exchange Contributions 

For the contributions to the matrix element of the third-order change to the reduced interaction from 
exchange processes, shown in Figs. 4-10, 4-11, and 4-12, we use 



\K) 



(130) 



and the change of variables 



p 4 = {xpf,xp 1± +r±), 

P5 = ([1 -x]pf, [1 -x]pi± -fj_), 



(131) 



which leads to the intermediate free masses 



M 2 L 



and 



x(l — x) 



Ml = ^ 



(1-x) 



(f± - q±) 2 

x-y 



(132) 



Pe = P2, 

P7 = {[x -y]pt, [x -y]pi± + [r± - q±\), 
Ps = P5- 



(133) 



Then 



+ 
+ 
+ 



1 9 

1Z AT f 1 

E (0203 1 6V^ |^>. = — J d 4 r ± dx-^—^— 



V) 



Ml 



KL 



Ml 



FK 



e 2A'- 4 A FK M|. 



1 

i 

^FK 



Ml 



Mf - 2Ml 



A FK Ml + A KL Ml 



^>-\A FK Ml+A KL Ml) 



FK 
1 



Ml 



FL 



„2A'- 4 A rr M? 



A 



KL 



Ml - 2M| 



2A'- 4 (A FJC Mi+A KI/ M^) 



[A' - A] 



} 



(134) 



This five-dimensional integral also cannot be done analytically. From Fig. 4, 

15 12 12 

E (^203 ^ (3) \<h)i = E fofal 6V(3) Mi = E <^2 3 | SV^ |0! 



i=13 



i=10 



2^3 



= 10 



(135) 



7.4 The Running Coupling 

We now have the complete matrix element of the third-order change to the reduced interaction in terms 
of numerical functions of q±, the transverse momentum of particle 2 in the center-of-mass frame, and y, 
the fraction of the total longitudinal momentum carried by particle 2. We first use the matrix element 
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to compute the scale dependence of g A . According Eq. (109), we need the p2± — P3± — limit of 
(</>2</>3| SV^ 0i). For both p2± and ps± to be zero, q± must be zero, in which case Mp is also zero. 
In this limit, all the integrals in (0203| SV^ 0i) can be done analytically, although care must be used 
because there are delicate cancellations of divergences among the different parts of each integrand. Note 
that the q± = limit of (0203 1 SV^ |0i) has no dependence on p^ or p^ . The result of applying Eq. 



(109) is 



B\,2 = C 3 = 



512^ l0g 



which from Eq. (W3) implies 



.9 A 



3 A' 2 
^ + 512^ 3 - l0g A^ 



A 2 

A 2 ' 



(136) 



(137) 



This equation tells us how the coupling is related at two different scales. Since A' > A, Eq. (137) shows 
that g A > g A , . This means the coupling grows as we reduce the amount by which free masses can change. 
This shows that the theory is asymptotically free, as expected. 

In conventional covariant perturbation theory with the minimal subtraction (MS) renormalization 
scheme, one obtains, for the scale dependence of the coupling [17], 

3 



,9m = V + 



5127T 



a' 2 



0(gl 



(138) 



where fj, and p! are two different renormalization points. We see that our coupling changes with our 
cutoff in the same manner that the MS coupling changes with the renormalization point. This is what 
one expects due to the scheme independence of the one-loop beta function. However, since there is no 
simple connection between our scheme and the MS scheme, we expect there to be no simple relation 
between the respective couplings at higher orders. 



7.5 Numerical Results for the Three-Point Matrix Element 

We now proceed to calculate the matrix element of the cutoff-dependent part of the third-order reduced 



interaction. From Eq. (11C) and the fact P2_l = P3± = implies that q± = 0, 



(02031^(^10!) = 



(020 3 |^ (3) |0l)-(020 3 |^ (3) |0l) 



<?1^0j A tor 



(139) 



where we use the fact that the q± = limit of (0203| SV^ |0i) is independent of p\ and p^ . Since 
the Hamiltonian is invariant under transverse rotations and q± is the only transverse vector on which 



(0 2 03 1 6VW |0i) depends, Eq. jl39| ) indicates that (0 2 3 | SV^ \<j> x ) is a function of q]_ . By inspection, we 
can see that (0203| SV^ |0i) is the product of N and some dimensionless quantity, and is the difference 
of a function of A and the same function with A — > A'. Also, by our assumption of approximate transverse 



locality, the RHS of Eq. (13£) can be expanded in powers of q±. Thus 



2 3 |5F (3) |0i) - (0 2 3 |^ (3) |0l) 



OO 



^V/*(y)9f (a 



-2i 



i=l 



A 



(140) 



where the fiS are unknown dimensionless functions of y. This implies 
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<f> 3 \V^(A)\fa) = 



(0 2 3 | 5V<® |0i) - (0 2 3 | SV& |0i) 



A terms 



i=l 



20 3 |5U (3) |0l}-(0203|^ (3) |0l} 



(141) 



Using this equation and our result for (02^3 SV^ \<f>i), it is straightforward to write (02^3 ^cd(^) I0i) 
as a numerical function of q±, y, and A. 

By inspection, (0203 |0i) is symmetric under y — > 1 — y and each (0 2 03| SV^ 3 ' |0i)j for i < 9 

is a function only of Mp, the free mass of the final state. It is not clear from inspection if the rest of 
{02 03 1 SV^ 3 ' \<pi) is a function only of M F and there is no reason that it should be. This means that the 
matrix element of the Hamiltonian is not necessarily a function only of Mp. 

Our result for the three-point matrix element of the Hamiltonian through third order, neglecting any 
cutoff-independent contribution to the third-order reduced interaction, is 

-A- 4 Mi " 



^,23^+^(02031^ (A) |0O 



(142) 



To get an idea of the size of the noncanonical contribution to the matrix element, we can compare it to 
the canonical contribution. To do this, we need to choose a value for the coupling. We would like to 
choose a large coupling so that we can get a pessimistic estimate of whether or not the expansion for the 



Hamiltonian is converging. When the coupling is large, the second term in Eq. (137) is nearly as large 
as the first. The value of log(A' 2 /A 2 ) is our choice, but the natural value is 1, because then the range 
of scales over which off-diagonal matrix elements of the Hamiltonian are being removed is comparable to 
the range of scales that remain. Thus we estimate that a large coupling is = 5127r 3 /3. 

Using this coupling, Fig. 5 compares the canonical and noncanonical contributions to the matrix 
element as a function of the free mass of the final state. The plot is the result of a numerical computation 
using the VEGAS Monte Carlo algorithm for multidimensional integration [18]. The contributions to 
the matrix element are plotted in units of the matrix element of the unregulated canonical interaction, 
and we consider only the situation in which the two final-state particles share the total longitudinal 
momentum equally. The solid curve represents the canonical contribution and the diamonds represent 
the noncanonical contribution. The statistical error bars from the Monte Carlo integration are too small 
to be visible. 

Fig. 5 shows that for the case we consider, the noncanonical corrections to the matrix element of 
the Hamiltonian are small compared to the canonical part, even when the perturbative expansion of 
the coupling is breaking down. However, this does not necessarily imply that corrections to physical 
quantities from the noncanonical part of this matrix element would be small. Investigation of this would 
require a fifth-order calculation in this theory; so we feel that this is not worth investigating until we 
consider QCD. 
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Figure 5: A comparison of the canonical and noncanonical contributions to (0203 1 M 2 {h) |0i) as a 
function of the free mass of the final state. The contributions to the matrix element are plotted in units 
of the matrix element of the unregulated canonical interaction, and we consider only the situation in 
which the two final-state particles share the total longitudinal momentum equally. The coupling used 
in the plot is g 2 — 512n 3 /3. The solid curve represents the canonical contribution and the diamonds 
represent the noncanonical contribution. The statistical error bars from the Monte Carlo integration are 
too small to be visible. 
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Figure 6: The noncanonical part of (02<fe| M 2 {h) \(j)i) as a function of the fraction of the total longitudinal 
momentum carried by particle 2 for fixed free mass of the final state. It is plotted in units of the matrix 
element of the unregulated canonical interaction for three different values of the free mass of the final 
state. The coupling used in the plot is g 2 = 5127r 3 /3. The statistical error bars from the Monte Carlo 
integration are too small to be visible. 



For the same coupling, Fig. 6 shows the noncanonical part of the matrix element as a function of the 
fraction of the total longitudinal momentum carried by particle 2 for fixed free mass of the final state. It 
is plotted in units of the matrix element of the unregulated canonical interaction for three different values 
of the free mass of the final state. This plot shows that the noncanonical part of the matrix element is a 
function of y for fixed Mp, and thus is not a function only of Mp. 



8 Conclusions 

We have outlined three conditions under which hadron states will rapidly converge in an expansion in 
free-particle Fock-space sectors in a HLFQCD approach. First, the diagonal matrix elements of the 
Hamiltonian in the free-particle Fock-space basis must be dominated by the free part of the Hamiltonian. 
Second, the off-diagonal matrix elements of the Hamiltonian in this basis must quickly decrease as the 
difference of the free masses of the states increases. Third, the mass of a free state must quickly increase 
as the number of particles in the state increases. We have argued that the Hamiltonian cannot meet these 
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conditions unless the vacuum is dominated by few-body free sectors; so we force the vacuum to be empty 
by requiring every particle to have a positive longitudinal momentum. We have argued that if the vacuum 
is empty and the Hamiltonian can be derived perturbatively, then the first and third conditions on the 
Hamiltonian should be automatically satisfied in HLFQCD, and the second condition can be enforced 
by suppressing the Hamiltonian's matrix elements between states that differ in free mass by more than 
a cutoff. 

The cutoff we use, like all regulators in HLFFT, violates a number of physical principles of the theory. 
This means that the Hamiltonian cannot be just the canonical Hamiltonian with masses and couplings 
redefined by renormalization. Instead, the Hamiltonian must be allowed to contain all operators that are 
consistent with the unviolated physical principles of the theory We have shown that if we require the 
Hamiltonian to produce cutoff-independent physical quantities and we require it to respect the unviolated 
physical principles of the theory, then the matrix elements of the Hamiltonian are uniquely determined 
in terms of the fundamental parameters of the theory. 

To force the Hamiltonian to produce cutoff-independent physical quantities, we require it to be unitar- 
ily equivalent to itself at a slightly larger cutoff. To make the Hamiltonian respect the unviolated physical 
principles of the theory, we place a number of additional constraints on it. First, the Hamiltonian must 
conserve momentum and be invariant under boosts and transverse rotations. Second, it should respect 
approximate cluster decomposition. Third, it should be approximately local in the transverse directions. 
Fourth, it should reproduce the correct free field theory in the nonintcracting limit. Fifth, it can depend 
only on the fundamental parameters of the theory and the cutoff. Sixth, the Hamiltonian must produce 
the correct second-order scattering amplitudes. The seventh and final restriction we place on the Hamil- 
tonian is that it must lead to perturbative scattering amplitudes with the correct dependence on the 
coupling. 

If the physical principles that we use are consistent with the physical principles violated by the cutoff, 
such as gauge invariance and Lorentz invariance, then the matrix elements of the Hamiltonian that follow 
from our procedure must lead to physical quantities that respect all the physical principles of the theory 

Our key assumption is that in an asymptotically free theory the matrix elements of the Hamiltonian 
can be computed perturbatively, provided the Hamiltonian is defined at a scale where the couplings 
are sufficiently small. This assumption and the rapid convergence of hadron states in the expansion in 
free-particle Fock-space sectors will allow us to approximate physical quantities in HLFQCD with a finite 
number of degrees of freedom. 

As an illustration of our method for calculating the Hamiltonian, we have calculated some of its 
second- and third-order matrix elements in (j> 3 theory in six dimensions. We have shown that the second- 
order matrix elements naturally obey the restrictions we place on the Hamiltonian and that they remove 
the cutoff dependence from physical quantities. We have also shown how the coupling depends on the 
scale in our approach and how asymptotic freedom arises. Finally, we have demonstrated that the third- 
order noncanonical corrections to the 4>\ — > 02<fe matrix element of the Hamiltonian that arise from the 
cutoff-dependent part of the reduced interactions are small, even when the perturbative expansion of the 
coupling is beginning to break down. 

We intend to apply our method to QCD; however, direct extension of the method to light-front QCD 
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is complicated by the presence of quark masses. This is because quark masses are additional fundamental 
parameters on which the Hamiltonian's matrix elements can depend. This complicates the method for 
determining the matrix elements, and thus the extension of our method to include particle masses requires 
further research. 

On the other hand, if quark masses are unimportant, then the only fundamental parameter in QCD 
is the quark-gluon coupling. It is then quite straightforward to apply our method to determine the 
Hamiltonian. There are two main complications that are absent in scalar theory. The first is algebraic 
complexity due to the presence of spin, color, and multiple canonical vertices. The second is the need for 
an additional cutoff to regulate divergent contributions to interactions involving gluons with vanishing 
longitudinal momenta. This cutoff can be implemented on the canonical interactions at the very beginning 
of the calculation of the Hamiltonian and can be taken to its limit at the end of a calculation of a physical 
quantity. This second cutoff does not complicate our method for determining the Hamiltonian, but having 
to take the cutoff to its limit significantly complicates the process of diagonalizing the Hamiltonian. One 
of the requirements that we have placed on the method presented here is that it will lead to physical 
quantities that are finite and well-behaved in this limit. A demonstration of this quality is part of the 
natural next step, application of our methods to a nonperturbative relativistic HLFQCD calculation. 
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APPENDIX A: Particle-Number Dependence of Free-Hamiltonian 
Eigenvalues 

The free Hamiltonian satisfies the Fock-space eigenvalue equation 

M) \F) = Mp \F) , (143) 

where \F) is any free-particle Fock-space state. We define Np to be the number of particles in \F). The 
eigenvalue Mp depends on Np: 



X 
i=l 



where Xi is the fraction of the total longitudinal momentum of \F) that is carried by particle i and f\± is 
the particle's transverse momentum in |f )'s center-of-mass frame. Momentum conservation implies that 

N F 

X>i = l. (145) 



If Np is increased, the x^s get smaller to maintain Eq. ( |145| ). Let us assume that 

1 



in order to maintain Eq. (145). Let us also assume that r?j_ is non-negligible and approximately inde- 
pendent of i. Then 

M f ~J2n f =N f , (147) 

i=l 



If the assumption in Eq. (146) does not hold, then Mp increases even faster with Nf. [j^] 



This argument fails if many particles have negligible ccntcr-of-mass transverse momenta. However, 
this is unlikely in a confining theory such as QCD because confinement suppresses any state containing 
particles with large transverse separation from the other particles in the state, and thus favors states 
containing particles with some non-negligible center-of-mass transverse momenta. 



11 Note that in equal-time field theory the free energy of a state increases only linearly with the number of particles, 
assuming they have non-negligible momenta. 
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APPENDIX B: Canonical Light-Front Massless Scalar Field Theory in Six 
Dimensions 



The purpose of this appendix is to state our conventions. A more complete discussion of canonical 
light-front scalar field theory can be found in Ref. [19]. 

With our conventions, any six-vector a is written in the form 



a — (a + , a , £Jj_)) 
where in terms of equal-time vector components, 



a ± =a°± a 5 



and 

4 4 

a± = ^2 a % ± ei = ^ a 4 e 4 , 

i=l i=l 

where e,; is the unit vector pointing in the i th direction. The inner product is 

a ■ b = -a + b~ H — a~b + — a± ■ b±, 
2 2 

and 



a_|_ = a± ■ a±. 



A spacetime coordinate is a six- vector, and according to Eq. (148), it is written 

x — (x + , x~ , x±). 



(148) 



(149) 



(150) 



(151) 



(152) 



(153) 



The time component is chosen to be x + . Here x~ is referred to as the longitudinal component, and x_\ 
contains the transverse components. 

The gradient operator is treated just like any other six-vector. Its components are 



and 



d ± = 2 



d 

d 

dx l . 



(154) 



(155) 



The canonical Lagrangian density for massless scalar field theory with a three-point interaction is 

1 f*A f>a 9 A3 



(156) 



2 r r 3! 

where the gradient operator is understood to act on the first spacetime-dependent function to its right, 
unless otherwise indicated by parentheses. The canonical Hamiltonian follows from C and the assumption 
that the field vanishes at spacetime infinity. It is 

H = I d A x\dx~ 



(157) 
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To this point, (j> has been regarded as a classical field. We work in the Schrodingcr representation, 
where operators are time-independent and states are time-dependent. Thus we quantize the field by 
defining it as a time- independent function of free-particle creation and annihilation operators: 



<j>{x ,x±) = J D x 



where 



and 



Di 



axe' 1 " 1 -* + a\e l - 



d 4 pi±dpf 



x+=0 



64tt 5 p+ 



Pi > o- 



(158) 



(159) 



(160) 



p + and p± are the conjugate momenta to x and x±; so they are referred to as the longitudinal and 
transverse momenta, respectively. The creation and annihilation operators follow the convention 



and have the commutation relations 



and 



a i = a(p~t,Pij_), 



[a,i,a\] = Si j 



[ai,a,j] = [at, at] = 0, 



(161) 



(162) 



(163) 



where 



and 



Sabc-,a>b>c>- - 647T 5 (p+ + p+ + p+ + • • -)<5 (5) ( Pa + Pb + Pc + Pa' - PV - Pc ) (164) 



SWfa-pj) = S(pf-p+)S^{p i± -p j± ). 



(165) 



Let V be the six-momentum operator and Ai be the invariant-mass operator. Since the momentum 
conjugate to x + is p~ , the Hamiltonian is identified as V~ , and it follows from 



V 2 = M 2 



that 



V~ = 



Pl + M 2 
V+ ' 



(166) 



(167) 



The Hamiltonian can be written as the sum of a free part and an interacting part: 

p-=H = h + v, (168) 

where 

(169) 



h= f D 1 fya\a 1 
J Pi 
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and 



v = (2tt) 5 . 9 / D X D 2 D S [ala ia2 5^ (p 3 - Pl - p 2 ) + 44<*i& 6) + P3 - Pi)] ■ (170) 

The process of normal ordering h and v produces infinite constants that have no physical significance 
and are dropped. The Hamiltonian also contains operators that have nonzero matrix elements only if 
a particle can have a longitudinal momentum of zero. These particle are dropped from the theory and 
thus the associated operators have no effect and are dropped. We plan to extend our method at some 
point in the future by replacing the effects of the dropped particles with interactions in the Hamiltonian. 
Without these particles, the vacuum is empty. 
The eigenstates of h are 



\4>i4>2 ■ ■ ■ 4>n) = a\a\ ■ ■ ■ al |0) 
for any integer n > 0. The associated eigenvalue equation is 

n 

h\</>i<fa---<l>n) = ^2p7 10102 •■•<; 

where 



i=l 



jT 2 

Pi± 



(171) 



(172) 



(173) 



and the sum is zero if n = 0. 



The nonintcracting limit of Eq. (167) is 



h = 



Ml 



where Mf is the free invariant-mass operator. It has the eigenvalue equation 



M) |0102 • • • 0n) = M 2 ■ ■ ■ 4> n ) 



where 



M 



Pi 



Pi 



(174) 



(175) 



(176) 



and P is the total momentum of the state. 

Finally, in terms of the free states, the completeness relation is 

1 = |0> (0| +/*£>! |0i) (0i | + i J D X D 2 |<^ 2 ) (0x0a | + 



(177) 
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APPENDIX C: Proof of Unitarity of U(A, A') 

[/(A, A'), as d 
T(A)f/(A, A') 



The purpose of this appendix is to prove that U(A, A'), as defined by 

dU(A,A') 



d(A-*) 



and 



f/(A,A) = l, 



(178) 
(179) 



is unitary as long as T(A) is anti-Hcrmitian. Assume that T is anti-Hermitian. Multiply Eq. (178) on 
the left by W(A, A'): 

^(A.AO ^y = C/t(A,A')T(A)C/(A,AO 

= - [f/ t (A,A')T(A)J7(A,A')] t 



1 1 



This is the same as 



^[f/t(A,A')£/(A,A')] =0, 



d(A"*) 

which implies that f/^(A, A')f/(A, A') is independent of A. Since Eq. ( |179| ) implies 

C/ t (A',A')L/(A',A') = 1, 

we conclude 

C/ t (A,A')C/(A,A') = 1. 



Now multiply Eq. (g_78|) on the right by U'(A, A'): 

^^V ( A , A = T(A)U(A, A')U^ (A, A'), 

and take the adjoint of this: 

U(A,A')^^1 = -C/(A,A')C/t(A,A')T(A). 



Adding Eqs. ( ]T84| ) and (|185|) gives 

(l7(A,A')l7 t ( A . A/ )) = [T(A),C/(A,A')C/ t (A,A')] 



d(A-*) 



FromEq. (179), 



U(A',A')U\A',A') = 1. 



(180) 
(181) 
(182) 
(183) 
(184) 
(185) 
(186) 
(187) 



U(A, A')U^(A, A') is a function of A that is uniquely determined by Eqs. ( |186| ) and ( |187| ). Therefore, 
since the statement 

tf(A,A')tft(A J A / ) = l (188) 

is a solution to Eqs. (186) and (187), it is a true statement. Since U(A,A') satisfies Eqs. (183) and ( |188| ), 
it is unitary. 
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APPENDIX D: Derivation of Reduced Interactions in Terms of Lower-Order 
Reduced Interactions 



In Section 4, we derived a constraint on the 0((f) reduced interaction for r > 1: 



(189) 



s=2 



Since we already know the first-order reduced interaction [see Eq. (pq)1, we wish to use this equation to 
compute the 0(g[) reduced interaction for r > 2, in terms of the lower-order reduced interactions. 

D.l The Cutoff-Dependent Part 

We begin by computing the cutoff-dependent part. Recall that momentum conservation implies that 
any matrix element of V(A) can be written as an expansion in unique products of momentum-conserving 
delta functions. This means that an arbitrary matrix element of Eq. ( |189|) can be expanded in products 
of delta functions: 

^(F|yW(A)|7)«-^(F|^(A')|/) W 



E 



(F\ SV^ \I) - ]T B r _ s>a (F\ V^~ S \A) \I) 



s=2 



(190) 



where the (i) superscripts denote that we are considering the i th product of delta functions that can 
occur in a delta function expansion of (F\ V^^A) |7). This equation is equivalent to a set of equations, 
one for each possible product of delta functions: 



(F\ 8VW \I) - J2 B r-s,s (F\ V^ S \A) \I) 



N 



F^({ Pn },A), 



(191) 



(192) 



W 

(F|^(A)|7)«-(F|yW(A')|/) (i) 
Cluster decomposition implies that we can write 

(F\V&(A)\I) {i) = 

where 5j is the j th delta function in the i th product of delta functions (it includes the longitudinal 
momentum factor), is the number of delta functions in the i th product, and F^({p n },A) is a 

function of the cutoff and the momenta of the particles in the matrix element, but does not contain delta 
functions that fix momenta. We define N p&It to be the number of particles in state |7) plus the number 
of particles in state \F), and n = 1, 2, 3, ... , N palt . Here p n is the momentum of particle n. We define 
N^ t to be the number of interacting particles in state |7) plus the number of interacting particles in state 
\F) for the i th product of delta functions. In order for the Hamiltonian to have the dimensions (mass) 2 , 
Fqj) must have the dimensions (mass) 6_2JV int . Note that we are suppressing the dependence of the RHS 
of this equation on r. 

We have assumed that any matrix element of the Hamiltonian can be expanded in powers of transverse 
momenta, not including the momentum-conserving delta functions; so 
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6-2N (i l 



E z 

{m nt } 



Wpart 4 



ipt}) UU( p -f 



(193) 



where t = 1, 2, 3, 4 denotes a component of transverse momentum and m„t is a non-negative integer index 
associated with transverse momentum component t of particle n. The sum is over all values of each of 
the rrint's, subject to the constraint that 



6 - 2JV« - $> nt ^ 0, 



(194) 



which is necessary to avoid terms in the momentum expansion that are cutoff-independent. The 2;| m "*^'s 
are the coefficients for the momentum expansion. They depend on the m nt 's and are functions of the 
longitudinal momenta of the particles. 



Since the RHS of Eq. (191) has the same product of delta functions as the LHS, we can write 



(F\ SV^ B rs,s ( F \ V(7 '~ S) ( A ) l J > 



G^({p n },A,A% (195) 



where G^ has dimensions (mass) 6-2 ^'^ , and by inspection of the LHS of Eq. ( |191| ) and Eq. ( |192| ) is a 
function of the momenta of the particles, A, and A'. Substitution of Eqs. ( |l95[) and ( |l92|) into Eq. ( |l9l[) 
yields 

Fgl({p n },A)-F%l({ Pn },A') = A, A'), (196) 

where the momenta in this equation are constrained by the delta function conditions. 

Since the LHS of Eq. (196) is the difference of a function of A and the same function with A — > A', 
G^ must be as well. Since the LHS of Eq. (jl9e| can be expanded in powers of transverse momenta, C?W 
must have the form 



G« ({*>„}, A, A') = £ ^ m " l} ({p+} 



{m„ t } 



A 



6-2JV. (i , 



n 



A 



-A 



I6-2N. 



n 



A' 



(197) 



where the sum is restricted by Eq. (194). 

Substituting Eqs. (|l97|) and (|l9g| ) into Eq. (|l9^) , we find 



A l> 



A 



/6-2JV. 



{m„t} 

E * 

{m, lt } 



n . / 



{"»„*} ( r + 



{m nt } 



A 



6-2AT. (i t ) 



n 



A' 



A 



A' 



n 

n,t 



A' 



(198) 
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Matching powers of transverse momenta on both sides of this equation gives 



J Y6-2Af i ( ^-S„, t m, lt _ ^/6-2N^-H n ,t m„t 
A 6-2JVW-S„, t m„ t _ A /6-2JV«-E„, t m n4 



(199) 



The factor in brackets cannot be zero because A ^ A' and Eq. (194) holds. Thus Eq. (199) implies 



r {m„ t } 



r{m nt } 



Then Eqs. (11931 ), ( |197|) , and (g00|) imply 



F^({p n },A) = G«({M,A,A') 



(200) 



(201) 



A terms 

where "A terms" means that is to be expanded in powers of transverse momenta and only the terms 



in the expansion that depend on A are to be kept. From Eqs. (192), (195), and (201), 



(F|^(A)|/)« = 



(F\ | J) - J2 B rs, s (F\ V^(A) \I) 



s=2 



(202) 



- A tern 



where it is understood that the momentum-conserving delta functions are ignored for the purposes of 
transverse-momentum expansions. Since a matrix element is the sum of the contributions to it from 
different products of delta functions, both sides of this equation can be summed over i to obtain 



(F\V^(A)\I) 



(F\ | J) - Br- S , s (F\ V^ S \A) \I) 



(203) 



A terms 



This equation tells us how to calculate the cutoff-dependent part of the O(g^) reduced interaction in 
terms of lower-order contributions. 



D.2 The Cutoff-Independent Part 

To complete the solution, we need to specify how to compute the cutoff- independent part. It is useful 
to first consider which contributions to V^ r \A) can be cutoff-independent. 

A matrix element of the cutoff-independent part of V^(A) can be expanded in products of delta 
functions and in powers of transverse momenta just as was done for the cutoff-dependent part. Thus we 
can write 

(F\V&\I)= y £{F\V&>\I}®, (204) 

i 

where 

F$({Pn}) (205) 
and 

■ATpart 4 * £ >. TTlnt 

= a 6 - 2 ^ y: ^(m) n n ( ¥ ) . ( 2oe ) 

{m„ t } n=l t=l ^ ' 



(F\V&\I) 



.(«) 



n4 

3=1 
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where the sum is over all values of each m nt , subject to the constraint 

6 - 27V.M _ J2 m nt = 0. (207) 

n,t 



Eq. (207) ensures that all the terms in the expansion of are cutoff- independent 



Eq. (207) places constraints on the possible cutoff-independent contributions to the reduced interac- 
tion. Any contribution to a matrix element of V^ r )(A) has an N^ t > 2, but Eq. (207) can only hold if 



-W/nt < 3. Suppose that = 2. In this case, F^l({p n }) can be written as a sum of terms, where each 
term corresponds to a distinct pair of interacting particles and depends on their momenta and the total 
longitudinal momentum [see Eq. ([32]) for an example]: 

rn 

= J2 F ci m \sm,P + ), (208) 

in 

where s m and s' m are the momenta for the initial and final particles in the m interacting pair, and 
where we have used the fact that for = 2, momentum conservation implies s m = s' m . According to 



Eq. ( |207| ), if N^. = 2, is cutoff-independent only if *£ n t m nt = 2. This means that either F^{ m ^ 
is quadratic in s m ± or it is zero. However, the matrix elements of the Hamiltonian are boost-invariant, 
as is the delta- function product in Eq. (205). This means that FqI must be boost-invariant, but it 



cannot be if Fq^ is quadratic in s m ±. Thus Fq^ = 0, FqI = 0, and there are no cutoff-independent 
contributions to the reduced interaction with two interacting particles. 

Note that contributions to the Hamiltonian with two interacting particles are self-energies, and they 
change the particle dispersion relation. If they change the dispersion relation such that the coefficients 
of the free relation become modified by interactions, then this can be viewed as renormalization of the 
field operators, i.e. wave function renormalization. This effect is absent unless either F^j or F^ can 
be quadratic in transverse momenta for N^ t = 2. We have just shown that this is not possible for i^n 



and according to Eq. ( |I94| ), Fq^ cannot be quadratic in transverse momenta for N^. — 2; so there is no 
wave function renormalization at any order in g A in our approach. 

The only other possibility for cutoff-independent contributions to the reduced interaction is from three 
interacting particles. Recall that in Section 3 we placed a restriction on the Hamiltonian that stated that 
each perturbative scattering amplitude that it produces must depend only on odd powers of the coupling 
if the number of particles changes by an odd number and only on even powers of the coupling if the 
number of particles changes by an even number. This implies that the O(g^) reduced interaction can 
change particle number by any even number s satisfying < s < r if r is even, or any odd number s 
satisfying I < s < r if r is odd. (We neglect to present the proof of this because it is tedious .P^]) This 

(r) 

means that since Vqj must have three interacting particles, it can be nonzero only if r is odd. 



According to Eq. (207), if N£ t — 3, then Fqj must be independent of all transverse momenta. Thus, 
if r is even, then = 0, and if r is odd, then is the part of V^ r \A) with three interacting 

particles and no transverse momentum dependence. Since VW is cutoff- independent, it should satisfy 
these conditions and inspection of Eq. (|3^ ) shows that it does. 

12 We assume that the frcc-statc matrix elements of the Hamiltonian can be chosen to be real. 
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(VI 

Before proceeding with the calculation of V^j , we need to deduce a bit more about the relationship 
of g A to g ' . According to Eq. (^lj) and the surrounding discussion, this relationship is determined by 
the matrix element (</>2<fe| SV \ <f>i), which can be expanded in powers of g , : 



b 2 fo\5V |0i> = Y^al {<h<h\sv {t) \<k) 



(209) 



t=3 



Recall that 6V& is built from products of V^(A')'s. This implies that SV^ can change particle number 
by 1 only if t is odd, and thus Eq. ( |209] ) implies that the coupling runs at odd orders; i.e. C s is zero if s 
is even [see Eq. (^3[)1. 

To calculate the cutoff-independent part of the 0(g r x ) reduced interaction, consider Eq. (191) with 
r -> r + 2: 



(F|^(A)|7)W-(F|T/^(A')|/) ( 



r+l 



{F\SV( r +V\I)-J2 B r+2-^ ( F \ V(7 



'(A)|/> 



6 = 2 



(i) 



(210) 



In order for l/( r )(A) to have a cutoff-independent part, r must be odd. We know V' 1 ' ; so assume r > 3. To 
calculate the cutoff- independent part of V <r - ) (A) , we only have to consider the case N.® t = 3. According 
to Eqs. p^ ), p3| ), and fll94] ), if A^ = 3, then for any t, 

^l^(A)|/>^., = 0, (211) 



where "Ext. fc^ — > 0" means the limit in which the transverse momenta in the external states are taken 



to zero. Then the zero-external-transverse-momentum limit of Eq. (210) is 



= 



(F\5V (r+2) | J) 



(0 

Ext. k^^O 



r+l 

B r +2-s,s 



(F\ V^ +2 - s) \I) W 



(212) 



s=2 



where the cutoff-dependent part of the matrix element in the sum is annihilated by the limit and the 
cutoff-independent part is untouched. If we take the first term in the sum on the RHS and move it to 
the LHS, we get 

r+l 



B r ,2 (F\ V£> \I) W 



(F\ SV {r+2) \I) 



Ext. k^^O 



^ Br+2-s.s {F\ V C j 



(r+2-s) | A (t) 



I) 



(213) 



We can sum over i, keeping in mind that we are restricted to the case of three interacting particles. Then 
we find 

r+l 



(F\V&\I) = 



1 



[F\ SV^+V \I) 



(r+2-s) 
CI 



Ext. k i ^0 s=3 



(214) 

3 int. part. 



where "3 int. part." means the products of momentum-conserving delta functions that appear on the 
RHS are to be examined and only the contributions from three interacting particles are to be kept. 

This equation tells us how to compute the cutoff- independent part of the O(g^) reduced interaction 
for odd r > 3 in terms of the cutoff-dependent part and lower-order reduced interactions. It is an 
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integral equation^] because Vqj is nested with lower-order reduced interactions in integrals in SV^ r+2 \ 
A cursory examination of the definition of SV^ r+2 ^ may lead one to believe that Eq. ( 214 ) is useless 
because it requires one to know I^ r+1 -'(A'). However, since r is odd, y( r+1 '(A') = V^ r D +1) (A), and the 
dependence of SV ( - r+2 ^ on V^d +1) (A') can be replaced with further dependence on Vq^ and V^q(A') (and 
lower-order reduced interactions) using Eq. (203) with r — ► r + 1. 



13 It is very difficult to prove that integral equations of this type have a unique solution; so we simply assume it is true in 
this case. 
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